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We perform realistic first-principles calculations of iron chalcogenides and ruthenate based mate- 
rials to identify experimental signatures of Hund's coupling induced correlations in these systems. 
We find that FeTe and K I Fe2- l/ Se2 display unusual orbital dependent fractional powerlaw behavior 
in their quasiparticle self energy and optical conductivity, a phenomena first identified in SrRuOa. 
Strong incoherence in the paramagnetic state of these materials results in electronic states hidden 
to angle-resolved photoemission spectroscopy which reemerge at low temperatures. We identify the 
effective low energy Hamiltonian describing these systems and show that these anomalies are not 
controlled by the proximity to a quantum critical point but result from coexistence of fast quantum 
mechanical orbital fluctuations and slow spin fluctuations. 



INTRODUCTION 

The study of the Hund's coupling effects in solids has a 
long history. Van der Marel and SawatskypQ pointed out 
that unlike the Hubbard U which is strongly screened, 
the atomic Hund's Jh, persists essentially unrenormal- 
ized in the solids and increases the splittings between the 
lower and the upper Hubbard bands for a half filled shell, 
while decreases it away from half filling. The Hund's 
term was also shown to have important consequences on 
the low energy physics of quasiparticles, when a tran- 
sition metal impurity is screened in metallic host. The 
Hund's coupling was shown to dramatically reduce the 
value of the Kondo temperature 2 , 3J. Recent interest in 
this problem arose from dynamical mean field theory [1] 
(DMFT) investigations of the recently discovered iron 
pnictide superconductors. It was proposed that in these 
materials strong correlations arise from the Hund's rule 
coupling Jff[5j rather than from the Hubbard U, result- 
ing in large mass enhancements. These calculations^ 
showed that for a reasonable value of the Hubbard U, 
the mass enhancement due to interactions is very small 
when Jjj=0, whereas it is exponentially enhanced by 
the Hund's rule coupling. Optical spectroscopy studies 
have shown that in both iron pnictides and chalcogenides 
the optical masses are many times larger than the band 
masses [SHI]. The trend in mass enhancements is well 
accounted for by DMFT combined with density func- 
tional theory (DFT+DMFT) calculations [10]. Since the 
strength of correlations in these solids is almost entirely 
due to the Hund's coupling, these materials are dubbed 
Hund's metals)]!]]. The role of Hund's coupling in iron 
pnictides and chalcogenides has been addressed from dif- 
ferent perspectives in the literature. [TUl - f2l)] Powerlaw be- 
havior in quasiparticle self energy of model Hamiltonians 
with Hund's coupling was discovered in Ref. [37] and re- 
lated to observations in ruthenates|28|l25]. Many anoma- 
lous properties of ruthenates [30] and other 4d compounds 
were shown to be governed by Hund's physics. [3"T] 

While at low energies and low temperatures Hund's 
metals are describable by Fermi liquid theory, the phys- 
ical properties in their incoherent regime are anomalous 
and surprising. In the iron pnictides and chalcogenides 



there is a strong tendency towards orbital differentia- 
tion [TU] , and the large mass enhancement can occur even 
though no clear Hubbard band exist in the one particle 
spectra of these Hund's metals |32 a . 

In this article, we use first principles methods and 
model Hamiltonians to search for experimental signa- 
tures of Hund's physics in iron chalcogenides and ruthen- 
ates which are the subject of current intensive experimen- 
tal studies. We show how the incoherence in iron chalco- 
genides above the Neel temperature can blur portions of 
the Fermi surface rendering them dark to photoemission 
spectroscopy. We show that fractional powerlaw behav- 
ior in optical conductivity that received significant atten- 
tion in the ruthenates also takes place in the FeTe sys- 
tem, deepening the analogies between these systems. We 
compare the powerlaw exponents in optical conductivity 
extracted from first principles DFT+DMFT calculations 
with experiments in a broad class of materials, and elu- 
cidate the control parameters that govern this behavior. 
The fractional powerlaw behavior is characteristic of an 
intermediate regime where the orbital degrees of freedom 
are quenched but the spin degrees of freedom are not. 
This physics is most pronounced at the special valence of 
one unit of charge away from the half-filling. 

RESULTS 
DFT+DMFT results 

We first show results of our realistic DFT+DMFT 
calculation |33l I34j for three compounds currently un- 
der extensive investigations: ruthenates Sr2Ru04, and 
iron chalcogenides FeTe and K a; Fe2_ a Se2. In all the 
DFT+DMFT calculations, we use the same Coulomb 
repulsion [7=5.0 eV and Hund's J#=0.80 eV, as de- 
termined from ab initio in Ref. 1321 and used in our 
previous work [TO] 112] . Notice that we include all 
the electronic states in a large energy window (-10 eV 
to 10 eV) as opposed to 3- or 5-band model calcula- 
tions, thus the Hubburd U appropriate for our calcu- 
lations is much less screened compared to those stud- 
ies. The electronic charge is computed self-consistently 
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FIG. 1: (color online) Fractional powerlaw in (a) theoret- 
ical self-energy and (b) experimental optical conductivity in 
iron chalcogenides and ruthenates. Experimental data are 
taken from ref. [32, for FeTe .9i and FeTeo.7Seo.3, ref. [33] 
for Fei.oeTeo.ssSo. 14, ref. [28] for Sr2Ru04, and ref. [29] for 
SrRu0 3 . 



on DFT+DMFT electronic density. The quantum impu- 
rity problem is solved by the continuous time quantum 
Monte Carlo (CTQMC) method [351 [36], using Slater 
form of the Coulomb repulsion in its fully rotational in- 
variant form. We use the experimentally determined lat- 
tice structures, including the internal positions of the 
atoms for Sr 2 Ru0 4 [37], FeTe [35], K x Fe 2 - v Se 2 [35], and 
FeSe @D]. 

We show in Figflja) the imaginary part of the 
self-energy of the t2g orbitals in Sr2Ru04, FeTe and 
KzFe2-ySe2 on the imaginary axis, plotted in login, scale. 
In the intermediate energy range from a low energy cut- 
off ^0.1 cV, below which the materials gain coherence, 
to roughly Hund's Jjf^0.8 eV, the imaginary part of the 
self-energy clearly shows a fractional powerlaw behavior, 
i.e., 7m£(iw„) cx — u>". For the normal Fermi liquid, this 
exponent is unity, and at finite temperature correlated 
materials have additional constant scattering rate. The 



fractional powerlaws are however very uncommon. 

From the quantum chemistry perspective, both iron 
chalchogenides and ruthenates share a common theme: 
they contain correlated electrons with the d valence of 
one unit charge away from half-filling. In iron pnic- 
tides/chalcogenides, the Fe-ion is surrounded by a tetra- 
hedron of pnictogen/chalcogen, and the resulting crystal 
field splittings are very small compared to Fe-pnictogen 
hybridization [S] hence all 5 Fe 3d orbitals are active. 
Their average occupancy is close to d 6 , one unit of charge 
away from the half-filled d 5 . For the ruthenates, the coor- 
dination of the Ru is octahedral, and the oxygen ligands 
induce a large t2g-eg splitting, with only the t2g orbitals 
active, containing approximately 4 electrons in three t2g 
orbitals, one electron more than the half filled shell. 

The values of the apparent powerlaw exponents differ 
from material to material and deviate even for different 
orbitals of the same material, which is connected to the 
orbital occupancy. As shown in Fig. [l|a), the xz/yz or- 
bitals of Sr2Ru04 show an exponent of 0.5, while the 
more correlated xy orbital, which is closer to half-filling, 
show a smaller exponent of 0.42 . In iron pnictides and 
chalcogenides, the average occupancy per orbital is even 
closer to half- filling (only 1/5 away, as opposed to 1/3 in 
ruthenates). As we will show below by a model study, 
one expects stronger electronic correlations in this case 
and a smaller powerlaw exponent. Indeed, the xz/yz 
orbitals in FeTe show exponent of ~ 0.36 whereas for 
the more strongly correlated xy orbital, the exponent is 
only « 0.24. Iron vacancies in the Ka;Fe2-ySe2 makes the 
compound even more correlated than FeTe, and the pow- 
erlaw exponent is further reduced to 0.27 for the xz/yz 
orbital, and only 0.07 for the xy orbital. 

The powerlaw behavior of the self-energy manifests it- 
self in optical conductivity studies. In a simplified treat- 
ment, the optical conductivity can be approximated by 
<t(w) oc Re(l/(u + »E"(w) + E'(w) - E'(w = 0))) gT]. In 
Fig. [TJb) , we present experimental data on FeTen qi |4"2"] . 
FeTeo.7Seo.3li2, Fei.oeTeo.ssSo. 14 03], Sr 2 Ru0 4 [2E], and 
SrRuCK[29]. As can be seen in Fig. [ljb), the optical con- 
ductivity in these materials can be roughly approximated 
by o\{ui) ~ Buj~ a in about the same energy range as in 
the theoretical self-energy. The experimental exponents 
obtained from optical conductivity are very similar to 
the theoretical exponents for the self-energy, as expected 
from the simplified relation between optical conductivity 
and self-energy. 

Hunds metals have a very low temperature scale, 
called the coherence temperature, below which a Fermi 
liquid-like coherence regime is reached. This phenom- 
ena has been discussed in other contexts such as heavy 
fermions|44| 45] and transition metal oxides [46] and can 
be fruitfully probed by photoemission spectroscopy. At 
finite temperatures, some electronic states can be very 
incoherent and coherence in different electronic states is 
usually not reached simultaneously. Due to the strong 
orbital differentiation discussed in Ref. [TO], the t 2g or- 
bitals have lower coherence temperature than the e g or- 
bitals in the iron-based superconductors. Within the t 2g 
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shell, the xy orbital has the lowest coherence temper- 
ature. In Fig. [2] we show the gradual evolution of the 
Fe t 2g orbitals from very incoherent state at high tem- 
perature to partially coherent state at lower temperature 
in paramagnetic (PM) state of FeTe. We display the 
momentum/orbital resolved density of electronic states 
at temperatures of 387 K, 232 K, 116 K and 58 K. 
For comparison, we also show momentum resolved den- 
sity of state of PM FeSe at 116 K, where all electronic 
states are quite coherent. The buildup of coherence in 
orbitally resolved spectra of FeTe is seen as a gradual 
buildup of the quasiparticle peak from a broad hump at 
elevated temperature to sharper peak at lower tempera- 
ture in Figs.[2je-f). In momentum space, the coherence is 
achieved more unevenly. While some bands can be iden- 
tified at 116K, and become pretty sharp at 58K, other 
bands are barely noticablc even at 58K. In particular, 
the band of primarily xy character circled by blue el- 
lipse, has enormous scattering rate at 58K and should be 
hard to be detected by ARPES. The missing Fermi sur- 
face is drawn in Fig. [2^ as large red pocket centered at 
r point, which is very incoherent above Tjv, hence miss- 
ing in the photoemission of the paramagnetic FeTe, in 
strong contrast to paramagnetic FeSe (Fig. [2k, see also 
Ref. [47]). Our calculation shows that K-intercalated 
FeSe (K x Fe 2 - y Se2) is even more correlated than FeTe, 
has smaller powerlaw exponents, and lower coherence 
temperature than FeTe. This is in agreement with recent 
angle-resolved photoemission spectroscopy (ARPES) ex- 
periments on A x Fe2- v Se2 compounds (A=K, Rb, Cs) 
where orbital dependent incoherence-coherence crossover 
was observed by Yi, Shen and collaborators |48j . 



Low energy Hamiltonian 

To gain some understanding of the Hund's physics in 
these systems, we derive below a low energy Hamiltonian 
of the three band Hubbard model, the simplest model 
which shows powerlaw behavior of the self-energy. The 
starting Hubbard Hamiltonian is H = Ht+Hjj, with the 
hopping term H t = £i 3 v, a ,& C/ioa/jto and Coulomb 

term ify = | £ i(Ti0bcd U l a > b > c > ^/L/L'/i^'iW Here 
a,b,c,d are orbital (site) indices, and a stands 

for the spin. The hopping term is taken to be locally 
SU(6) symmetric (no crystal fields), while the Coulomb 
interaction is set to U[a,b, c, d] = U5 a dSbc + J5 a c5bd, 
which reduces the symmetry to SU(3)xSU(2). Within 
DMFT, this model maps to a SU(3)xSU(2) impurity 
Hamiltonian. To understand why the Hund's rule 
coupling has such a dramatic effect on the physical 
properties, we first perform Schrieffer- Wolff trans- 
formation (for derivation see appendix) to obtain 
Kondo-like Hamiltonian, of the form H^°J ldo — 
Hq + Hi + H2 + H3 , with the potential scattering term 
Ho = •/ P £ a(T V'L( )^w( )> the spin-spin Kondo part 
Hi = ^lE^EwVi^OK^VwCO), the orbital- 
Kondo part, H 2 = J 2 £ Q T a E ol)(7 ^ a(T (0)A Q al ,^(0), 
and the coupled spin-orbital part H$ = 
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FIG. 2: (color online) Incoherence-coherence crossover in 
FeTe. A(k, u)) of along the path Y^X^M^Y^Z^ 
R -> A -> Z for FeTe at (a) 232 K, (b) 116 K, and (c) 58 K 
and for (d) FeSe at 116 K in the PM states, (e) and (f) A(oj) 
for Fe 3d xy and xz/yz orbital at 387 K, 232 K, 116 K and 
58 K in PM FeTe. (g) and (h) color-coded Fermi surface in 
the F plane for PM FeTe and FeSe, respectively. Red, green 
and blue color correspond to xy, xz and yz orbital character, 
respectively. Due to the incoherent nature of the xy orbital 
above Tn, the outer hole pocket around Y is not easy to be 
detected in ARPES experiment. 



h Eq/3 T a SP J2aba ^a«m a ab<J^«^b*'(0). Here 

S a = Eacrcr' faa^oo'faa' an d T@ — Eabcr facr^abfba' 

are spin, and SU(3) orbital operators acting on the 
impurity site, tp(0) are field operators of the conduction 
electrons coupled to the impurity, while a^ a , and A" & 
are Pauli matrices and the Gell-Mann 3x3 matrices of 
the SU(3) group, respectively. 

Notice that in our picture the same electrons carry 
both orbital and spin degrees of freedom, in contrast to 
the point of view of Ref. I49L which emphasizes the spin 
and orbital degrees of freedom being carried by different 
type of electrons, i.e., tig the spin, and eg the orbital. 

While the form of the low energy impurity model is 
dictated by symmetry considerations, the exchange cou- 
plings Ji , J2 , J3 depend crucially on the impurity valence 
and Hund's coupling Jh- For the half-filled shell and 
large Jh, only the spin-spin term J\ survives, and a well 
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known reduction of the J\ Kondo coupling for a factor 
of (21 + 1) was derived in Refs. |2j [3] compared to a corre- 
sponding one band model. Consequently, a huge reduc- 
tion of the Kondo temperature for a factor of (21 + l) 2 
in the exponent was derived in Ref. 3. This regime is 
relevant for half-filled d 5 shell realized in the Hund's in- 
sulators LaMnPO [50] . 

For the above Hund's metals, the relevant valence of 
transition metal ion is one unit of charge away from half- 
filling. When Jh is negligible, the Hamiltonian is SU(6) 
symmetric, and all three Kondo couplings Ji, J2, J3 are 
positive (antiferromagnetic). For the valence n imp = 2 
(or riimp — 4), their numerical values are J\ = Jo/3, 
•h = Jo/4, and J 3 = J /2, where J = V 2 /(2U + e f ) 
(or Jo = V 2 /(3U + £/)) is a positive number, which de- 
pends on the corresponding Anderson impurity model 
parameters, i.e., hybridization V and impurity level et. 
The ground state is a Fermi liquid, because antiferro- 
magnetic couplings between conduction electrons and 
impurity degrees of freedom ensure complete quenching 
of both the orbital and spin moments. On the other 
hand, when Jh is large, the spin-spin Kondo coupling 
Ji changes sign to ferromagnetic, while the orbit J2 
and spin-orbit J3 remain positive. In the three band 
SU(2) x SU(3) model and large Jh, their numerical values 
are J\ = — Jo/9, J2 = Jo/3, and J3 = Jo/3, where Jo = 
V 2 /(2U-2J H +e f ) > (or J = V 2 /(3U+ J H + £/) > 0) 
for valence n imp = 2 (or n,; mp = 4). This change of sign 
is due to the orbital blocking mechanism [TU] , that allows 
only those virtual charge excitations which go through or- 
bital singlet intermediate state (see appendix). We note 
that for valence d 6 in iron pnictides/chalcogenides, the 
correct low energy Hamiltonian has three terms, but not 
just spin-spin term, as proposed earlier |51j . It is how- 
ever the spin-spin J\ term that changes sign in the limit 
of large Jh, and impedes quenching of the spin degrees 
of freedom (termed spin freezing in Ref. 27). This sub- 
stantially reduces coherence temperature, however, the 
J3 term, which couples spin and orbital, is positive and 
gives rise to the Fermi liquid state at very low tempera- 
ture. 



Model Hamiltonian calculations 

To demonstrate the above picture, we numerically 
solve a simplified three band model with the nearest 
neighbor diagonal hopping t aa — 0.4, and the next 
nearest neighbor hopping t' aa — 0.4 and i! a <g = 0.2, 
which gives a total bandwidth of the tight banding model 
W 1=3 3.5. We take U = 6 and large Hund's coupling 
J/f=2 and Jh — 1 for powerlaw to extend over a larger 
frequency range. 

In Fig. [3^,, we show the imaginary part of the quasi- 
particle self energies for Jh=2. At the intermediate tem- 
perature T=0.01, the self energies of n,j=2.00 and 2.26 
display finite values at zero frequency by extrapolation, 
which suggests incoherent properties at this temperature. 
However, at a lower temperature T=0. 00125, the corre- 
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FIG. 3: (color online) Quasiparticle self energies for a three 
band model with Jh = 2.0 at different filling rid=1.75, 2.00, 
and 2.26. (a) The self energies at two temperatures T=0.01 
and 0.00125 show the incoherence-coherence crossover with 
decreasing temperature, (b) The self energies plotted in logio 
scale display fractional powerlaw behavior in the intermediate 
frequency range from £q to Ei as indicated by arrows. 

sponding self energies clearly display Fermi liquid behav- 
ior at low frequencies. Therefore there is an incoherence- 
coherence crossover with decreasing temperature. We 
determine the coherence temperature as the tempera- 
ture at which the renormalized scattering rate is equal 
to the temperature, i.e., ~zImY,(iu = i0 + ,T*) = k^T* , 
where 1/z = 1 — dImY,(iu = i0 + ,T*)/duj and ks is the 
Boltzmann constant. We show the coherent temperature 
in Fig. |4^ as a function of electron occupation rid- We 
reached 8-times lower temperature than previous stud- 
ies [27] to access the Fermi liquid state at filling far be- 
yond rid — 2, and map out the coherence incoherence 
crossover temperature T* . 

A good powerlaw fit to the self-energy, as shown in 
Fig. [3]d, is obtained only in a limited range of frequency 
between the low energy cutoff proportional to the Fermi 
liquid scale (eq in Fig. [3|d), and the high energy cutoff 
(el), which is always smaller than the Hund's coupling. 
The range of frequencies at which the powerlaw is valid 
(sq < Lo n < £*, as indicated by arrows in Fig. [SJd) is 
largest at valence rid = 2, where the exponent is close 
to 1/2, as previously reported in Ref. [27]. Much lower 
temperature reached in this work show that exponent a 
decreases monotonically with increasing rid (see Fig. [3 d 
and Figs. [5] and [6] in appendix), in contrast to Ref. [27 , 
hence stronger correlations approaching the Mott state 
at rid — 3 lead to smaller exponent at low temperature, 
in qualitative agreement with our realistic calculations 
for Fe and Ru compounds. Most importantly, there is 
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FIG. 4: (color online) (a) The coherent temperature; and (b) 
the spin and orbital susceptibility, as a functional of rid for a 
three band model with Jh = 2.0 (solid lines) and Jh = 1.0 
(dash lines). 



no signature of quantum phase transition to non-Fermi 
liquid spin- frozen state around valence rid — 2, where 
the powerlaw exponents are found, and the crossover 
temperature scale does not follow the powerlaw behavior 
T* ~ (x — x c ) z,y expected in a quantum critical scenario. 
Hence we can exclude the possibility that the exponents 
are due to the proximity to the quantum phase transition, 
as proposed in Ref . [27l |52] . 



DISCUSSIONS 



The powerlaw exponents are found in the temperature 
and frequency regime where the spin degrees of freedom 
are very slow and unquenched (spin susceptibility has 
Curie- Weiss form and large static value) while the orbital 
degrees of freedom are very fast and quenched (orbital 
susceptibility is Pauli-like, but enhanced) . (see Fig. [4jz> 
and Fig. [7] in appendix) This is a novel regime in which 
two degrees of freedom behave in different ways , one 
fluctuates very fast (positive J2), the other one fluctuates 
very slowly (negative Ji), and both coupled by a third an- 
tiferromagnetic coupling J3. This situation is similar to 
the intermediate phase of the extended Hubbard model 
study of Ref. |53j , where similar continuously varying ex- 
ponents were shown to exists in the metallic non-Fermi 
liquid phase in which there was a quenched spin-degree 
of freedom and an unquenched charge degree of freedom. 
Notice that at rid = 3 only the large spin state is pos- 
sible, hence the orbital degrees of freedom are gapped, 
and exponents disappear, while the effect is maximal one 
unit of charge away from half-filling, i.e., rid — 2. 



CONCLUSIONS 

In conclusions, we have shown in this paper that the 
Hund's rule coupling has a strong impact on the elec- 
tronic states in valence of one unit of charge away from 
half-filling. The strongly correlated state in such mate- 
rials can have very low coherence temperature, and the 
self-energy and optical conductivity show fractional pow- 
erlaw at intermediate energy. We have derived the ef- 
fective low energy Hamiltonian describing these systems 
and identified a negative Kondo coupling in the spin-spin 
part of the corresponding low energy Kondo model. We 
have shown that these anomalies are not controlled by the 
proximity to a quantum critical point but result from co- 
existence of fast quantum mechanical orbital fluctuations 
and slow spin fluctuations. This is relevant for ruthenates 
and iron chalcogenides, as well as many other materials 
with similar valence and sizable Hund's coupling. 
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APPENDIX 
Schrieffer- Wolff transformation 

We start our discussion with the three band Hubbard 
model H — H t + Hu, with the hopping term 



Ht — J]] tabfiaafjber 
ija,a : b 



(1) 



and the Coulomb repulsion term 

H u = \ E ^M,M/L/L'iW'jw (2) 



Here index a runs over the three orbitals, i, j over lattice 
sites, and a over spin. The hopping term is taken to be 
locally SU(6) symmetric, while the Coulomb interaction 
is set to U[a, b, c, d] = US a dSbc + JSacdbd, which reduces 
the symmetry to SU(3)xSU(2). 

Within the Dynamical Mean Field Theory, this model 
maps to the SU(3)xSU(2) impurity Hamiltonian of the 
form H imp — Hbath + Hh y b + Hi oca i 



Hbath — E £ kaV'kaCT'/'ka<T 



kao" 



H h yb = V" Vka1pk aa faa + h.C. 



local 



katr 



£ ffaafaa + H L 



(3) 



(4) 



(5) 



G 



To gain further insights into the low energy degrees of 
freedom of this Hamiltonian, we perform the Schrieffer- 
Wolff transformation, which takes the form 



H, 



eff 



(6) 



where P n is the projector to the impurity ground state 
multiplet, and AE is the energy cost for the charge exci- 
tation from the ground state multiplet to the nil lowest 
energy multiplet states, and is always negative. 

Here we will limit our discussion to the case of a ground 
state valence rif = 2 and virtual charge excitations into 
valence n/ = 3. A direct way of evaluating this effec- 
tive Hamiltonian is to perform exact diagonalization of 
Hiocai , then expressing matrix elements of f a<T operators 
in terms of atomic eigenstates {F\ a ) mim2 — (mi\f\ cr \m2) : 
and finally evaluating all terms which appear in the sum 



H eff - 



kk' abaa' 



E ( F aa)m im2 (F^,) m2 m 3 \m 1 )(m 3 \. 



(7) 



Here m 2 runs over the ground state multiplet at valence 
n + 1, while mi and m 3 run over ground state multiplet 
at valence n. Finally we need to express the impurity 
degrees of freedom in terms of the impurity operators 
such as is the total spin S and the orbital isospin operator 
T. 

This tedious derivation can be circumvented by a trick. 
We Fourier transform the bath operators V 2 ip^ a (0) — 
^ka^kac an( l introduce the combined spin-orbit index 
i = (acr). We can then rewrite the effective Hamiltonian 
for charge excitations from valence n to valence n + 1 as 



\ " 



V 



ijkl 



4(0)^(0)P„/ fc P„ +1 //P„ 6(i, k)8(j, 0(8) 



We next find a complete orthonormal basis in the space 
of spin and orbit degrees of freedom (Tr(/ t V 3 ^) = <5 Q( g), 
in which the completeness relation takes the form 



S(i,k)S(j,l) 



xi a *)iAn 



kl- 



(9) 



Here matrices I a form a complete basis for the 
SU(3)xSU(2) group. For the SU(2) and SU(3) group 
we use Pauli 2x2 matrices <r a , and Gell-Mann 3x3 
matrices A". In terms of these, the complete basis I a is 




Next we realize that even in the presence of an arbitrary 
projector, the local operators keep the same form of the 
expansion in terms of the electron field operator 



J2(l®° a )kl PnfkPn + lf}Pn OC -S° 
kl 

YJ(A Q (8 l) fc , P n f k P n+1 fJP n cx -T° 



ki 



^(A Q ® <jP) kl P„/ fe P n+1 //P n cx -T° 



(12) 
(13) 
S", (14) 



kl 



but the proportionality constants need to be determined 
by an explicit calculation. Notice that just like in the 
Wigner-Eckart theorem, we only need to consider one 
matrix element to determine proportionality constant, 
which greatly simplifies this derivation. 



Now we can recognize that the first term in Eq. ( 10 ) 
gives rise to potential scattering of the form 



(15) 



the second term in Eq. ([lO]) gives the spin-Kondo part 

h x = .h ]T s a v>L(oK^Vw(o), (i6) 

a aaa' 

the third gives orbital-Kondo part 

H 2 = J 2 Ta E ^WUAV^OV (17) 

a otafi 

and the last gives spin-orbit Kondo part 

H 3 = J 3 E T a ® S b ^ aa X a a ^a b a ^(0)^ 

a, b aa/3 

(18) 

The Kondo-couplings J\ , J 2 , J3 depend on the valence 
rif and type of the projector P. We first consider the 
SU(6) symmetric case, which is realized in the absence 
of Hund's rule coupling. In this case, the projector P 
is irrelevant, since all states at some valence have equal 
energy. The local operators are then simply given by 

E faa&aa'fla' = ~ ^S*" 
aaa' 

Ef \ Q ft _rpa 
Jaa^abJba' ~ ± 

aba 



The normalization factors come from the fact that 
Tr(cr a cr^) = 28(a,P) and Tr(A Q A' 3 ) = 2<5(a,/3). We can 
then simplify the low energy Hamiltonian by 

Heff = E P n.fkPn+l.f}Pn 4 (0) ^ (0011) 

ijkl,a 



(10) and the Kondo couplings become J\ — 2/6 Jo = Jo/3, 
■h = Jo/4: and J 3 = 2/4 J , where J = 2 J^ e/ > 0. 
Notice that all Kondo-couplings are positive (minus sign 
comes from AE and from the proportionality constant) 
and hence antiferromagnetic couplings ensures complete 
quenching of the spin and orbital moment. The ground 
state is thus Fermi liquid. 

In the limit of large Hund's coupling, the projector 
Pn+i projects to the subspace of high-spin states only, 
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which in the case of the three band model and n/ 
take the following form 



|1) 
|2) 



13) = 
|4) = |ttt) 



|||) 

-Jg<Utt> 



iti> + 1 m» 
tit) + 1 m» 



3. 

(19) 
(20) 

(21) 
(22) 



Projection to the ground state multiplet P n is achieved 
by projecting to the following states 



|5) 
|6) 
|7) 

|8) 



Ho) 
in) 

|12) 
|13) 



l|0|) 
|0 ||) 
I ||0) 

1 

71 

1 
i 

71 
I tot) 

lott) 
Itt o) 



(Itot) 
(|oit) + 
(I it o) - 



I- It 04)) 
|0t4» 

■in o)) 



(23) 
(24) 
(25) 

(26) 
(27) 

(28) 

(29) 
(30) 
(31) 



An explicit calculation can be used to determine pro- 
portionality constants 

^ ] Pnfaa^cra'Pn+lfaryiPn = 7^5 
aaa' 

^ PrLfaaKib°aa'P'n+lft a :Pn = ~ ~T a 
abcra' 

We can finally determine the Kondo couplings in the limit 
of large Hund's coupling. Their value is J\ — —2/3 * 
1/6J = -Jo/9, J 2 = 4/3 * 1/4J = Jo/3 and J 3 = 
4/3 * 1/4J = Jo/3. Here J = 2U ^j H +e f > °- 

The crucial result of this calculation is that the spin- 
spin Kondo coupling J\ changes sign when Hund's 
coupling is strong. This comes from the fact 
that the spin operator in the projected subspace 
Ea.*'Pnfa.<T^,Pn + ifl,Pn = §S« has very different 
expansion in terms of electron field operator than in non- 
projected case Y,aaa> J 'av^a* /<L' = - 2 S a - The origin of 
this sign change is in the orbital blocking mechanism, 
which ensures that the intermediate state at nt — 3 is 
a high spin-state (in this case 5 = 1) but is orbitally a 
singlet state, such as states [l9p2l Orbital blocking is 
a restriction in the Hilbert space imposed by the large 
Hund's rule coupling. It modifies the Kondo couplings 
away from their SU(N) symmetric values (Ji = Jo/3, 
J2 = Jo/4, J3 = Jo/2). This blocking results in different 
Kondo couplings in different valences. For the half-filled 
shell (relevant for Mn 2+ ) it results in J 2 = 0, J3 = 0, 



and a strong reduction of the value of J\ , first recognized 
by Schrieffer [2J. For the valence of one unit of charge 
away from half-filling (relevant for Fe 2+ and Ru 4+ ), or- 
bital blocking results in the sign reversal of Ji. 



Results for the three band Hubbard Model 




FIG. 5: (color online) Quasiparticle self energy for 
Jtr=2.0. The imaginary part of the quasiparticle self en- 
ergy in logio-logio scale as a function of electron occupation 
rid for U =6.0 and Jir=2.0. Note the data are shifted along the 
y-axis for better illustration. The linear dispersion of the self 
energy in the plot indicates the powerlaw behavior exists in 
the intermediate frequency region as indicated by the arrows. 

Using our numerical methods of quantum Monte Carlo, 
we can not obtain high precision real axis self-energy, 
however, we can infer its analytic properties from imag- 
inary axis analogs. The fractional exponent in scat- 
tering rate on the real axis (/mE(w) cx u a ) leads to 
the same powerlaw on imaginary axis for the imaginary 
part (IrnS(iui n ) oc — w"). The real part, on the other 
hand, shows the powerlaw only when scattering rate is 
very asymmetric around zero frequency. For example, 
E"(w > 0) = A\x\ a and £"(cj < 0) = B\x\ a , the real part 
on imaginary axis is ReE(iui n ) cx (A — B) J Q A x\x\ a / (x 2 + 
ojfydx (where A is the upper cutoff for the powerlaw), 
and does not show powerlaw in the symmetric A = B 
case. Our calculation shows that the real part does not 
show very clear powerlaw on imaginary axis, hence we 
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FIG. 6: (color online) Quasiparticle self energy for 
Jff=1.0. The imaginary part of the quasiparticle self en- 
ergy in logio-logio scale as a function of electron occupation 
rid for U =6.0 and Jff =1.0. Note the data are shifted along the 
y-axis for better illustration. The linear dispersion of the self 
energy in the plot indicates the powerlaw behavior exists in 
the intermediate frequency region as indicated by the arrows. 

infer that the scattering rate is quite symmetric at low 
frequency on real axis. 

Figure [5] and [6] show the imaginary part of the quasi- 
particle self energy in logio-logio scale as a function of 
electron occupation rid for [7=6.0 and J#=2.0 and 1.0, 
respectively. The linear dispersion of the self energy in 
the plots indicates the powerlaw behavior exists in the 
intermediate frequency region as indicated by the ar- 
rows. For both values of the powerlaw exponent 
a decreases monotonically with increasing rid towards 
half filling, i.e., n^=3. The upper energy cutoff e\ drops 
rapidly for rid >2.0, suggesting the powerlaw behavior 
is vanishing quickly when rid goes away from 2 to half 
filling. Compared to Jff=2.0, the powerlaw behavior for 
J/f =1.0 is valid in a smaller frequency region and the cor- 
responding powerlaw exponent is larger, suggesting the 
important role of Hund's coupling in giving rise to the 
powerlaw behavior. Therefore the powerlaw behavior is 
most visible at electron occupation one unit of charge 
away from half-filling, in this case, rid = 2 

Figure [7] shows the local spin, orbital and charge sus- 
ceptibility at zero frequency as a function of temperature 
for n d =1.75, 2.00, 2.26 and [7=6.0, J H =2.0. The spin 



FIG. 7: (color online) Local spin, orbital, and charge 

susceptibility The local spin, orbital and charge susceptibil- 
ity at zero frequency as a function of temperature for n,j=1.75, 
2.00, 2.26 and [7=6.0, Jh=2.0. The spin susceptibility has 
large static values and takes the Curie- Weiss form while the 
orbital susceptibility is Pauli-like and enhanced at interme- 
diate temperature and around n,j=2.0. Note the charge sus- 
ceptibility is two orders of magnitude smaller than the orbital 
susceptibility thus doesn't play an important role. 

susceptibility has large static values and takes the Curie- 
Weiss form while the orbital susceptibility is Pauli-like 
and enhanced at intermediate temperature and around 
71^=2.0. Note the charge susceptibility is two orders of 
magnitude smaller than the orbital susceptibility thus 
doesn't play an important role. 
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